Analytical expression for the convolution of a Fano line profile with a Gaussian 
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Asymmetric Fano line profiles are frequently encountered, e.g., in the photoionization spectra of 
atoms and ions. For the fitting of spectral line profiles to experimental spectra the line profiles have 
to be convoluted with the experimental window function. The latter is often taken to be a Gaussian. 
It is shown that the convolution can be represented by a rather simple analytic expression involving 
the complex error function. 
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I. INTRODUCTION 

Asymmetric line profiles are frequently encountered, 
e.g., in atomic pfiotoionization due to a quantum me- 
chanical interference between resonant and a nonresonant 
ionization pathways. According to the quantum mechan- 
ical analysis of Fano [1] the photoionization cross section 
as a function of photon energy E in the vicinity of the 
resonance energy £"1.68 can be represented as 
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where ctq and ui are slowly varying functions, 
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is a reduced energy, A^ is the natural (Lorentzian) line 
width, and q is the asymmetry parameter. The latter 
can be calculated theoretically from the transition ma- 
trix elements of the interfering resonant and nonresonant 
ionization channels. 

As an example Fig. 1 shows the photoionization cross 
section of Be-like B"*" ions in the photon energy range 
26.5-31.2 eV which was measured at a synchrotron ra- 
diation source [2]. The resonances are associated with 
dipole allowed double photoexcitation of the ls^2s^ ^S 
ground state to ls^2pns ^P and ls^2pnd^P states and 
subsequent autoionization. The principal quantum num- 
ber n depends on the photon energy. For n — )> cxd both 
Rydberg series converge towards the series limit at about 
31.15 eV. 

Clearly the individual resonance line shapes are 
strongly asymmetric. For the extraction of the resonance 
parameters cti, E'res, ^l and q from the experimental 
data the resonance lines can be fitted by a Fano profile. 
In such fits the experimental photon energy distribution 
(window function) has to be taken into account. In many 
cases the experimental window function can be repre- 
sented as a Gaussian where the Gaussian full width at 
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FIG. 1. Rydberg series of resonances in the photoionization 
cross section of Be-like B+{ls'^ 2s^ ^S) ions [2]. The fuU line 
is a fit of Fano profiles convoluted with a gaussian to the 
experimental data points. The 2pns^P and 2pnd^P series 
have different asymmetry parameters (cf. Eq. 1) of ~ —0.4 
and ~ — 2, respectively. 



half maximum (FWHM) corresponds to the experimen- 
tal energy spread. Thus, for a fit to the experimental 
data the Fano profile has to be convoluted with a Gaus- 
sian. 

The situation is similar to emission spectroscopy of hot 
gases where Doppler broadening results in Voigt line pro- 
files, i.e., the convolution of a Lorentzian with a Gaus- 
sian. It is well known (see, e.g., [3-6]) that the Voigt 
profile can be calculated efficiently from the complex er- 
ror function [7] (also known as the Faddeeva function). 

It is much less known that also the convolution of a 
Fano profile with a Gaussian can be represented by the 
complex error function. A corresponding rater simple 
formula has already been given in the appendix of a pre- 
vious publication [8]. The formula allows for a fast and 
accurate evaluation of the convolution in peak fitting rou- 
tines. Here its derivation is presented. 

It is noted that a different, more complex formula has 
been published earlier without its derivation [9] . It seems, 
that this formula has not received much attention since 
even in more recent work the convolution of a Fano profile 



with a Gaussian has only been treated approximately or 
by numerical integration (see, e.g., [10, 11]). 

The present paper is organized as follows. In section 
II the calculation of the Voigt profile from the complex 
error function is reviewed. In section III some relevant 
properties of the complex error functions are elucidated. 
An expression for the convolution of the Fano profile with 
a Gaussian in terms of the complex error function is de- 
rived in section IV. A conclusive summary is given in 
section V. 



II. THE VOIGT PROFILE 

The convolution of a Lorcntzian line profile 
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with a Gaussian 
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yields the Voigt profile 
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The profiles in Eqs. 3 and 4 arc normalized such that 

L{E)dE = A 'And G{E)dE = l. (6) 

The widths A^ and Aq are the Lorentzian and Gaussian 
FWHM, respectively. With the definitions 
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Eq. 6 transforms into 
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where w{z) denotes the complex error function and z 
X + iy. 



III. SOME PROPERTIES OF THE COMPLEX 
ERROR FUNCTION 

The complex error function is defined for 5z ~ y > 

as [7, 12] 
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For later use we now calculate the integral 

dt. 
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For the calculation of the derivative on the right-hand 
side of this equation we exploit the identity Wriiz) = 
w{y/rjz) which follows from the substitution ix/rjt) — > t 
in Eq. 15. This yields 
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To this end we define 
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where in the last step we have used the identity [7, 12 
dw{z) 
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(8) Combining Eqs. 16 and 17 yields 
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From equation 18 the partial derivatives of w{z) with 
respect to x and y are easily calculated as 
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These are of use, e.g., in least-squares fitting routines 



IV. CONVOLUTION OF A FANO PROFILE 
WITH A GAUSSIAN 

For the purpose of peak fitting we define the Fano line 
profile somewhat differently as suggested by Eq. 1, i.e. 
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The term -1 inside the square brackets ensures that 
F{E) -^ ioT E ^ ±00. For q — > cx) the Fano profile 
as defined by Eq. 22 approaches a symmetric Lorentzian 
(Eq. 3), i.e.' F{E) -> L{E). With the definitions of i, x 
and y from Eqs. 7, 8, and 9, respectively, the convolution 
with a Gaussian as defined by Eq. 4 can be expressed as 
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In the equations above the definitions of 3fJ(w) (Eq. 12), 
3(m;) (Eq. 13), and I1 (Eq. 14) as well as Eq. 19 have been 
used. After gathering all terms one finally arrives at 
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It is easily seen that C{E) — > V^(i?) for g — > 00 as ex- 
pected. 



V. SUMMARY 

It has been demonstrated that the convolution of a 
Fano line profile (Eq. 22) with a gaussian (Eq. 4) can be 
evaluated analytically (Eq. 24) by using the complex er- 



I 

ror function wi^z) with z = x-\-%y and with x and y from 
Eqs. 8 and 9, respectively. Various fast and accurate al- 
gorithms for computing the complex error function have 
been described in the liteature [3, 4, 6, 13-17]. Their 
performances have been critically evaluated repeatedly 
[5, 6, 18, 19]. According to the findings of Zaghloul and 
Ali [6] their algorithm is the most accurate available to 
date. 

The line profile C{E) has been implemented by the 
author as a user-supplied fit function for the commercial 
software Origin [20] and successfully used in various 
contexts (see, e.g.. Fig. 1). The implementation is 
available from the author upon request. 
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